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Abstract 

The linearised time dependent coupled Gross-Pitaevskii equations describing 
the long wavelength excitations of Bose-condensed binary mixtures are solved, 
in the bulk and in harmonic traps for the case where only the binary phase is 
present. In the former case we obtain two zero-sound branches In the latter 
case the dispersion law also contains two branches whose dependence on the 
quantum numbers of the modes is the same as for a one-component conden- 
sate, but with different prefactors, depending on the ratios of the three s-wave 
scattering lengths of the two atomic species. In the general case where the 
binary phase in the trap coexists with one or both one-component phases, the 
mode spectrum depends on the geometry of the interphase boundaries due 
to the boundary conditions there. Measurements of the oscillation frequen- 
cies as in recent experiments with modulated traps would yield very detailed 
information on this system. 
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The recent dramatic advances in the ability to create Bose condensates have now reached 
the stage where Bose condensates of binary mixtures of two atomic species only distinguished 
in their angular momentum quantum numbers can be manufactured in a magnetic trap at 
effectively zero temperature The ground state of such binary Bose- condensed mixtures 
have recently been studied theoretically in the Thomas- Fermi approximation 0. To learn 
more about their properties it seems interesting to study their long wavelength excitations. 
Recent experiments for one-component condensates have been reported in |3|,§]], where col- 
lective excitations have been excited by changes in the trap potential. Theoretical work on 
collective excitations, again for a single- component condensate, has been presented in 
where the linearised Gross-Pitaevskii equation was solved numerically, and in ||, where an 
analytical solution was obtained in the long wavelength limit. Very good agreement between 
the theoretical predictions for the low lying mode frequencies [||,|| and the experimental re- 
sults |3|,|[] was obtained. As shown by Stringari H, in the Thomas Fermi approximation, 
the frequencies of the long wavelength excitations of a trapped Bose-Einstein condensate are 
independent of the microscopic properties of the condensate, in particular of the two-body 
scattering length, and only depend on the frequency u)o of the single-particle oscillation in 
the trap. 

In the present paper we shall consider the two coupled linearised Gross-Pitaevskii equa- 
tions at zero temperature describing the excitations of the binary phase of a mixture of two 
Bose condensates. Considering the spatially homogeneous system first we obtain the Bo- 
goliubov spectrum for this case with two branches of zero sound for k — > 0, where the sound 
velocities depend on the ratios of the three different scattering lengths in the problem. We 
also derive the coupled linearised hydrodynamic equations and their boundary conditions 
which govern the long wavelength excitations of Bose condensed binary mixtures in a trap 
coexisting with the one-component condensates. For the special case of a pure binary-phase 
condensate in a spherically symmetric trap we determine the mode-spectrum analytically 
in the long-wavelength and Thomas-Fermi approximations. We find that the dependence of 
the mode spectrum on the radial and angular momentum quantum numbers n and I is the 
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same as in the case of a single- component mixture, 



00+ - = Vb+,~ (2^ 2 + 2nZ + 3n + I 
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but the prefactor u; 2 + _ is no longer simply given by the trap frequency, Uq, but it has two 
solutions, corresponding to the two sound branches in the bulk, each of which depends on 
the two ratios of the two-particle scattering lengths. 

The coupled Gross-Pitaevskii equations describing a binary mixture of Bose-Einstein 
condensates with equal particle masses are given by 



Hx = -2^V 2 Vi + V>i (Wil 2 + W2I 2 ) + £Ma#i, 



hip 2 = -— vV 2 + ^2 (v 2 \H 2 + Vu\H 2 ) + u 2 (x)^ 2 . 
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We shall assume in this paper that the conditions 



(2a) 
(2b) 



V 1 >o,v 2 > 0, v x v 2 - v? 2 > 
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are satisfied. Here V* 



di, with i — 1,2, 12, are the interaction parameters, a, being 



the s-wave scattering lengths. Let us first consider a spatially homogeneous system, i.e. the 
case of vanishing trap potentials U\(x) = U 2 (x) = 0, where the solutions of (§) of the form 



ipi ~ e are given by 
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The Hi, with i = 1,2, are the chemical potentials of the two particle species to be 
determined by normalising \tpi\ 2 to the particle numbers N^. The stability regions of the 
three solutions have been discussed in D. In particular, the solution (2) for the binary 
phase is stable if its right-hand sides are both positive. 

Putting now 
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& = \f\W*¥+fre~ i{mt/h ~ (pi \ (5) 

and linearizing Eqs. (0) in pj and <pj we obtain for the excitations in the stability region of 
the binary condensate, 

<Pi = -t {ViPi + V 12 p 2 ) + - , TTiaT^ VVi, (6a) 
" 4m |^i | 2 

p 1 = --|^ 12) | 2 VV, (6b) 
m 

and similarly for <p 2 ,p2 simply by interchanging indices. The dispersion law for solutions 
(fi, pi proportional to e^i^-k-x) j g eas jiy obtained as to 2 = _{k 2 )k 2 with 

Ai,-(* 2 ) = ^ + <-, (7a) 



(2m) 
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For large k the two branches of the dispersion law approach that for the free particles like 

hk 2 mci _ , N 

w= + — 8 

The constant shift from the free-particle dispersion law would show up e.g. in light scattering 
as discussed for the case of a single component condensate in . For k — > 0, on the other 
hand, we obtain two zero-sound branches of collective excitations with sound velocities c +i _ 
where the densities of the two components oscillate out of phase and in phase respectively. 

Let us now turn to the case of a trapped binary condensate, Ui(x), t/ 2 (a;) ^ 0, which we 
consider in the Thomas-Fermi approximation. Then the time-dependent solutions (f|) still 
apply with the replacements pi —>■ pi — Ui(x)(see 0). Linearising Eqs. (Q) with the ansatz 
(||) and making a long-wavelength approximation by keeping only the lowest order spatial 
derivatives we obtain in the place of ([]) 

0i = -t (^iPi + ^i2P 2 ) , (9a) 

p 1 = --V-|# 2) (*)| 2 V<pi, (9b) 
m 
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and similarly for ip 2 , P2 by interchanging indices. Eliminating the phases ifii we obtain the 
two coupled equations for the long-wavelength density fluctuations, 

1 



p 2 = -V • \4 12) (x)\ 2 V(V 2 p 2 + V 12Pl ). (10b) 
m 

In the general case, where the [12]-phase coexists with the one-component phases k = [1], [2] 



h = -V • \^r>(x)\ 2 V(V lPl + V 12 p2), (10a) 
m 

m 



Eqs. (10) have to be solved simultaneously with the equations 

A = -V • |V>?Ww,a (ii) 

m 

and boundary conditions. These are dictated by number conservation and require the con- 
tinuity of pi and the normal derivative of tpi across the interphase boundary between the 
phases [12] and k on which ip^ ^ 0. In general, the spectrum of collective excitations of the 
total system will therefore depend on the details of the geometry of the interphase bound- 
aries, and solutions of the coupled equations can only be obtained numerically. Here we shall 
consider a special case where no interphase boundaries are present, for which an analytical 
solution can be given similar to the case of a one-component condensate considered in ref. 
||: Let us consider the case where Ui = ^mQfr 2 represent two concentric isotropic harmonic 
traps. It follows from the results in that in the special case where p 2 /pi = (Q 2 /Qi) 2 with 

,v 2 v 12 . n 2 v 2 v 12 . 

max(—,—) > -3 > min(—,—) (12a) 

V12 Vl S 4l V12 V\ 

no inter-phase boundaries occur in the trap and only the binary phase is present there. 
The condition for the chemical potentials requires that the particle numbers are chosen 
in the ratio N 2 /Ni = (ViOjj ~ ^i2^i)/(^2^i — ^12^2)' which can be satisfied mathemat- 
ically if ( p]2|) holds, but can be satisfied physically only if the ratios of the three scat- 
tering lengths is known. Under these conditions we can solve Eqs. ( |T0| ) by an ansatz 
Pi = e^ t T,k=oCik rl+2k Yim(9,(i>)- The coefficients C ik , i = 1, 2, k = 0, 1, • • • , n, are deter- 
mined by a recursion law. The condition Ci n+ \ = determines u 2 by Eq. ([]]) with 
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|^^a 2 (ai - ai 2 ) + fi2 a i( a 2 - ^12) 
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(f2 2 a 2 (ai + ai 2 ) - fi 2 ai(a 2 + ai 2 )) + 4a 2 2 (a 2 fi 2 - a^Vt^j (a^l - ai 2 fi 2 
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For ai 2 = the dispersion law (|i~3T) reduces to Stringari's result cj 2 + = Vt\ 2 for each 
particle species ||. The result (|T3|) is plotted in fig. la (upper branch) and fig. lb (lower 



branch). As long as the inequalities (12) are satisfied the squared frequencies (13) are both 
positive. At the borders of this region in parameter space the frequency vanishes and 
ujb + approaches Q± or f2 2 depending on which of the two borders has been reached. This 
can be seen happening in figs.la,b, where oJb + = Q 2 at the stability border defined by the 
vanishing of Ub- for the choice Vt\ = 2Qf made in these figures. 
For axially symmetric trap potentials, 



O 2 



cos 2 ^ 



(14) 



with fijL/fip = fij_ 2 /^±i special solutions of Eqs. ([U]) can still be obtained, as in the single 
component case ||. Solutions oc r l Yi m (9,(f)) can be found for 



m = ±/ : ujj 



+,-> 



(15) 



where 



and 



uj-. 



% li 



(16) 



m = ±(Z-l): cj 2 _ = (/ - + ai 



+,-> 
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where 



0J\\ 



b+,- 



n 2 ^n 2 



(18) 



For the four coupled m = modes with n = 0, I — 2 and n = 1, I = one obtains the four 
frequencies, 
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"l-i = "l + _ ( 2 + -Xl _ + - A /9Ai_ - 16Ai_ + 16 



2 A+ - + 2 

<-, 2 = (2 + - ~ 16A^,_ + 16) , (19a) 

with A+ _ = cj| + _/ojj_ + _. 

Generalisations of our results for mixtures with more than two components can be easily 
made. All that changes is the expressions for the prefactors io\ in the dispersion law (HI). 
E.g. for condensates of ternary mixtures in isotropic concentric potentials Ui(x) one would 
obtain three branches Eq. (P, whose squared frequencies w^ 123 are the roots of a cubic 
equation and depend on the ratios of all six scattering lengths ai, a 2 , a 3 , ai 2 , ai 3 , a 23 . 

Mathematically, all explicit results ([!]), (|I5D , (|T7j), (P^|) are the same as the results for the 
one-component condensate, if there the oscillation frequency Uq in the trap is replaced by the 
two frequencies flTBp. Physically, this difference is of interest, because these frequencies are 
determined dynamically by the coupling between the two components in the binary phase. 

Measuring the frequency of the low lying excitations of the binary condensate phase in 
coexistence with one or both single-component phases with basically the same experimental 
set-up as employed in Refs. HH seems very feasible. Because of the coupling of the coexisting 
phases by particle number conservation the mode spectrum is predicted to depend in general 
on the geometry of the interphase boundaries. To check our prediction for the simple mode 
spectrum in the pure binary phase the following requirements have to be met: The ratios 
of the scattering lengths must be known and inequality (|1^) must happen to be satisfied for 
the atomic species under consideration. Then in principle a pure binary phase condensate 
can be produced in concentric spherically or axially symmetric traps assuming that the 
ratios of the particle numbers can be chosen appropriately. To obtain concentric traps, their 
potentials must be made sufficiently stiff to suppress the offset of their centers due to gravity. 
Satisfying all these conditions one should be able to observe that the mode spectrum reduces 
to the simple form ([!]) with the scale factors (13). 
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FIGURES 

FIG. 1. The frequency scale factors u^/fif (figure l.a) and a^_/f2f (figure l.b) as a function 
of the ratio of the scattering lengths 012/01, for different values of 02/01, for f^/^i = 2. The 
chosen value of the ratio f^/^f corresponds to the experiment Q. 



10 



